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Abstract 

BES data on J/^ — ^{K^K^ir^) are presented. There is a strong peak due 
to ?7(1440)/i, which is fitted with a Breit-Wigner amphtude with s-dependent 
widths for decays to K*K, kK, rjirn and pp; n refers to the K-k S-wave. At 
a KKiT mass of ~ 2040 MeV, there is a second peak with width ~ 400 MeV; 

= 0~ is preferred over 1+ and 2" respectively by 5.2 and 6.8 standard 
deviations. It is a possible candidate for a 0~ ssg hybrid partner of 7r(1800). 

PACS numbers: 14.40.Cs, 12.39.Mk, 13.25.Jx, 13.40.Hq 
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Typeset using REVT^ 



There have been earher data from Mark III and DM2 0] for J / ip radiative decays to 
K^K^-K^ , as well as K^K'tc^. Recently, the BES group has published data on the latter 
channel 0. Here we present BES data on decays to K^K^-k^. These data have lower 
backgrounds than for K~^K~'tt^, because of the identification of Kg tt+tt^. Consequently, 
the partial wave analysis may be extended up to a KKtt mass of 2300 MeV, covering an 
interesting structure at ~ 2040 MeV. 

The Beijing Spectrometer(BES) has collected 7.8 x 10^ J/ip triggers, used here. Details 
of the detector are given in Ref. 0. We describe briefly those detector elements playing 
a crucial role in the present measurement. Tracking is provided by a 10 superlayer main 
drift chamber (MDC). Each superlayer contains four layers of sense wires measuring both 
the position and the ionization energy loss (dE/dx) of charged particles. The momentum 
resolution is cxp/P = 1.7%vl + P^, where P is the momentum of charged tracks in GeV/c. 
The resolution of the dE/dx measurement is ~ ±9%, providing good tt/K separation and 
proton identification for momenta up to 600 MeV/c. An array of 48 scintillation counters 
surrounding the MDC measures the time-of-flight (TOP) of charged tracks with a resolution 
of 330ps for hadrons. Outside the TOP system is an electromagnetic calorimeter made of 
lead sheets and streamer tubes and having a z positional resolution of 4 cm. The energy 
resolution scales as as/E = 22%/ V^, where E is the energy in GeV. Outside the shower 
counter is a solenoidal magnet producing a 0.4 Tesla magnetic field. 

Each candidate event is required to have four charged tracks. Each track must have 
a good helix fit in the polar angle range —0.8 < cosO < 0.8 and a transverse momentum 
> 60 MeV/c. A vertex is required within an interaction region ±30 cm longitudinally 
and 3 cm radially. A positive identification of just one is required using time of flight 
and/or dE/dx. Events are fitted kinematically to the 4C hypothesis J/ip —>■ 7(fC^7r^7r+7r^), 
requiring a confidence level > 5%. 

Backgrounds arise mainly from 11^X^71^71^71^ and K'^71^71^71^ . Those events giving a 
better fit to these channels are rejected. Next, we require | Umiss \ = \ Emiss — Pmiss \< 0.15 
GeV/c^, so as to reject the events with multi- photons or more or less than one charged 
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kaon; here, Emiss and Pmiss are, respectively, the missing energy and missing momentum of 
all charged particles. The momentum of the K^Ti^Ti^'n-^ system transverse to the photon 
Pl^ = 4 I Pmiss I ^ sin^(6'm^/2) < 0.005 {GeV/cf is required in order to remove the 
background J/ip ^ tt^K^tt^ii^h^] here Om^y is the angle between the missing momentum 
and the photon direction. Finally, are selected with a cut on the tt+tt^ invariant mass, 
I M^+^- - Mko \< 25 MeV. Fig. |l](a) shows the TT+TT invariant mass closest to the Kg 
mass before the K'^ are selected; a very strong signal Kg is seen. The number of surviving 
events is 1095 with 57 ± 5 non-Kg background under the Kg. For our final fit, we use 683 
events below a KKtc mass of 2.3 GeV. A constraint to the Kg vertex does not improve the 
signal/background ratio further, but loses some events. 

The effects of the various selection cuts on the data is simulated with a full Monte Carlo 
of the BES detector including the decay path of the K^; 250,000 Monte Carlo events are 
successfully fitted to J/ip ^ 7(ii'^ii"°7r^). All background reactions are similarly fitted to 
this channel. The estimated background is 29 ± 7%, mostly from J/ip ^ -k^ {K^ K^tt"^) , 
some from non-i^^^ events. It peaks at about 2.3 GeV, and follows phase space closely. We 
have included it in the amplitude analysis, but it has little effect, since all genuine signals 
have a characteristic dependence on either or both of production and decay angles. 

Fig. |l](b) shows the K'^K^Tt'^ mass spectrum; the dark histogram shows the estimated 
background in the analysis region. There is a conspicuous and somewhat asymmetric peak 
due to ?7(1440)/t, similar to the earlier data from Mark III, DM2 and BES. At high mass, 
there is a distinct peak at 2040 MeV. Fig. ^ shows Dalitz plots for three mass ranges: (a) 
1360-1560 MeV, (b) 1600-1750 MeV, and (c) 1800-2200 MeV; fits are shown in (d), (e) 
and (f). There is a conspicuous K*K decay mode in the first region of the ?7(1440). At 
higher masses, it disappears rapidly, and the mass projections shown in Fig. ^ are consistent 
with decays to k,K only; above 1560 MeV, there is no significant evidence for ao(980)fC, 
ao ^ KK. 

We have carried out a partial wave analysis using amplitudes constructed from Lorentz- 
invariant combinations of the 4- vectors and the photon polarization for J / ip initial states with 
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helicity ±1. Cross sections are summed over photon polarisations. The relative magnitudes 
and phases of the amplitudes are determined by a maximum likelihood fit. We include KKir 
states with quantum numbers 0^, 2" and 2+. There are two helicity amplitudes for 1+, 
three for 2~ and three for 2+. Because production is via an electromagnetic transition, the 
same phase is used for different helicity amplitudes to the same final state. Different phases 
are allowed for different decay channels, e.g K*K and kK, because of strong interaction 
effects due to rescattering. 

The analysis is discussed separately for the mass region of r7(1440) and the 2040 MeV 
peak. The 77(1440) has been fitted using a Breit-Wigner amplitude with s-dependent width: 

/ = ^ . (1) 

The numerator A is a complex coupling constant. The T{s) are taken to be proportional to 
the available phase space for each channel, evaluated numerically p. The 777177 phase space 
is taken from rja, the dominant channel, but ao(980)7r phase space is similar and both are 
slowly varying over this mass region. The magnitude of each F is adjusted iteratively so 
that cross sections integrated over the resonance agree with the branching ratios determined 
experimentally. The magnitude of F^o- has been obtained from BES data on radiative decays 
to ?77r7r . That for Fpp has been obtained by fitting BES data on radiative decays to 47r [0 , 
including in the fit 77(1440) and the broad 77(1800). Here //(ISOO) refers to the very broad 
0~ signal (F ~ 1 GeV) derived by Bugg and Zou from an analysis of several channels of 
J /ip radiative decay. Values of Tk*k and F^i^ are obtained from the present data. 

In the mass region of the 7^(1440), half the nK signal comes from the low mass tail of 
r;(1800) and its constructive interference with 7/(1440). That is, if 7/(1800) is removed from 
the fit, the kK width of 7/(1440) needs to be doubled. Removing the 7/(1800) has a significant, 
but not dramatic, effect on log likelihood, which changes by 4.8 for 2 extra parameters. 

The /i(1420) is also included in the amplitude analysis, and a small component due to 
/i(1285). Both optimise close to the masses and widths quoted by the Particle Data Group 
(PDG) [P], so we fix them at PDG values. The amplitude analysis distinguishes cleanly 
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between quantum numbers l"*" and 0~ for K*K decays. If the whole ?7(1440) signal is fitted 
with = 1+ (optimising its mass and width), log likelihood is worse by 11.4, a significant 
amount. (Our definition of log likelihood is such that it increases by 0.5 for a one standard 
deviation change in one parameter). 

In the earlier analysis of BES data on the K^K^tt^ final state Q], a fairly large amplitude 
was fitted for 7^(1800). The smaller background in present data and the wider mass range 
allow us to show that this component should in fact be rather small. Its effects on the 
7^(1440) may be replaced with some increase in the total width of that resonance and an 
increase in its width for decays to K*K. Present results for the fitted widths are shown in 
Fig. ^ and branching fractions in Table 1. The fit is compared with the KKtt mass spectrum 
by the histogram in Fig. 5. 

A free fit to the mass gives 1440 MeV. However, rjTnT data give a resonance mass of 
1405 ± 5 MeV, according to the summary by the Particle Data Group . The s-dependent 
width we use for 7/(1440) explains naturally a mass difference of 20 MeV between rjmr and 
KKtt data; the rapidly increasing phase space for K*K makes the KKn channel peak higher 
and explains also the asymmetric shape of the peak, which rises rapidly on the lower side of 
the peak and falls more slowly on the upper side. A small (~ 15 MeV) discrepancy remains 
between the peaks fitted to r/vrvr and KKtt. We adopt a compromise between fitting these 
data and i]tttt by using a mass of 1432 MeV, but the effect on other conclusions is negligible. 
Interferences between 7/(1440) and the broad ?7(1800) depend on their relative phases and 
can shift the peak in different data sets; so we do not regard this small discrepancy as a 
matter for concern. 

Around 1650 MeV, there is some indication for a narrow KKtt peak. However, fitting it 
requires an unreasonably narrow width ~ 30 MeV. An ss state at this mass has no obvious 
non-strange partners. If fitted, it is only a two standard deviation effect. Therefore we 
discard it as a statistical flutuation. Including it has neghgible effects on parameters fitted 
to r/(1440) and the peak at 2040 MeV. 

We now turn to the latter peak. It cannot be explained by the very broad r/(1800), 
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which has a completely different and much flatter shape, illustrated by the shaded area in 
Figure |^(b) below. We fit it with a simple Breit-Wigner amplitude of constant width. Its 
mass and width optimise at M = 2040 ± 50 MeV, T = 400 ± 90 MeV. We have tried fits to 
this peak with resonances having quantum numbers 0~, 1"^ and 2~; for standard qq states, 
one does not expect S"*" in kaonic channels until 2300 MeV. We find that log likelihood is 
better for 0~ than 1+ by 16.4. The latter has one additional parameter, so it is a poorer fit 
by 5.2 standard deviations. If a combination of 0~ and 1"*" amplitudes is used, log likelihood 
improves only by 0.6, and the fitted 1"^ component is very small: 4.4% of 0~ in cross section. 
These results are not sensitive to the //(ISOO) contribution: removing it, the distinction 
between quantum numbers 0~ and l"*" for the 2040 MeV peak remains at a log likelihood 
difference of 12.9. 

We have also tried adding or substituting 2^. Alone it gives a poor fit, worse in log 
likelihood than 0~ by 27.9. This demonstrates that 2~ and 0~ are well separated by their 
distinctively different angular distributions. If it is added freely to the fit, it improves log 
likelihood by 3.2 for three extra parameters; this cannot be considered significant. Fig. |^ 
shows magnitudes of components fitted in the amplitude analysis when the 2040 MeV peak 
is fitted as 0^. The slight differences between Figs. 5(b) and (d) is due to interferences of 
r/(1440) and r7(2040) with the broad r/(1800) in Fig. 5(b). 

Branching fractions for production and decay, including the dominant interferences, are 
given in Table 2. Values are integrated up to a KKtt mass of 2.3 GeV. Decays to K^Kstt^ 
have a branching ratio 1/3 of all KKtt decays. We correct all measured branching ratios 
by this factor 3, so as to quote branching fractions for all KKn charge states. The overall 
branching fraction, summed over all final states is (6.0 ± 0.4 ± 2.1) x 10^^. 

We now discuss possible interpretations for the 2040 MeV peak. Our data for J/ip 
radiative decays to r/vr+Tr^ [6] were fitted using an ?7(1760) with a width of 250 MeV and 
an 772(1840). The 7^(1760) is entirely distinct from 7/(1800), which has a much larger width. 
A possible interpretation is that it is the n = 3 qq state. Then the 7/(2040) observed here 
could be its ss partner. 
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However, the VES collaboration has identified a 7r(1800) [|T0|] with curious decay modes 
to /o(1300)7r, /o(980)7r and -K'o(1430)i^', but not pn. There has been speculation that this 
is an / = 1 hybrid |jTl|. The ?7(1760) would make a natural partner; its decays to rja and 
ao(980)7r are to be expected for a hybrid. It is natural to expect a corresponding ssg state 
decaying to kK in the KKtt channel roughly 200-250 MeV above the peak in rjini. In J/ if: 
radiative decays, the amplitude for production of qq states is suppressed by two powers of 
Os, required to couple intermediate gluons to quarks; at 2040 MeV, as — 0.41. Production 
of a hybrid will only be suppressed by one power of as in amplitude. We therefore examine 
the possible interpretation of 77(1 760) as a qqg hybrid. 

For a hybrid, the branching fraction expected in the KKn channel is half that for rjirir, 
since in J/ip decays intermediate gluons couple equally to uu, dd and ss. If fitted as 0~, 
the branching ratio for the 2040 MeV peak in KKtt is (2.1 ± 0.1 ± 0.7) x 10"^; this value is 
obtained after allowing for interferences with 7^(1800) and includes the error in the overall 
normalisation. It is to be compared with the branching ratio for ?7(1760) in rjinr of (1.8 ± 
0.75) X 10^'^ 1^. These values are consistent within the sizable errors with the expectation 
for hybrids. 

The magnitude of branching ratio we now fit to ?7(1800) kK is (0.58 ± 0.03 ± 0.20) x 
10^'^. Again, the error includes the overall normalisation uncertainty. It compares with 
(1.08 ± 0.45) X 10~^ fitted to titttt decays |0. Within the errors, these values are now 
consistent with fiavour-blind decays of a glueball. 

In summary, present data contain less background than earlier data on J/ip — >■ 
'^{K'^K'TT^) and allow a somewhat improved determination of the properties of 7^(1440). 
Its dominant decay mode is to K*K. This suggests it is the first radial excitation of 7^(958), 
probably mixed with the broad 7^(1800), in order to account for its strong production in J/ip 
radiative decays. We now find a small component of 77(1800) decaying to kK. 

We observe a peak at 2040 MeV which may be fitted with a 0~ resonance of width 400 
MeV. = 0^ is preferred over 1"*" and 2~ respectively by 5.2 and 6.8 standard deviations. 
Its branching fraction, when compared with the r]7C7c channel, would be consistent with 

8 



interpretation as a ssg hybrid. 

The BES group thanks the staff of IHEP for technical support in running the experiment. 
This work is supported in part by China Postdoctoral Science Foundation and National 
Natural Science Foundation of China under contract Nos. 19991480, 19825116 and 19605007; 
and by the Chinese Academy of Sciences under contract No. KJ 95T-03(IHEP). We also 
acknowledge financial support from the Royal Society for collaboration between Chinese and 
UK groups. 



9 



REFERENCES 

^ Data analyzed were taken prior to the participation of U.S. members of the BES Col- 
laboration. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7 

[s; 

[9 

[lo: 
[11 

[12 

[13 



J.Z. Bai et al, Phys. Rev. Lett. 65 (1990) 2507. 

J.-E. Augustin et al, Phys. Rev D46 (1992) 1951. 

J.Z. Bai et al., Phys. Lett. B440 (1998) 217. 

BES Collaboration, Nucl. Instr. Methods, A344 (1994) 319. 

D.V. Bugg, A.V. Sarantsev and B.S. Zou, Nucl. Phys. B471 (1996) 59, equn. (40). 

J.Z. Bai et al, Phys. Lett. B446 (1999) 356. 

J.Z. Bai et al., Phys. Lett. B472 (2000) 207. 

D.V.Bugg and B.S.Zou, Phys. Lett. B 396 (1997) 295. 

Particle Data Group, Euro. Phys. J. C3 (1998) 1. 

D. Amelin et al., Phys. Lett. B356 (1995) 595. 

F.E. Close and P.R. Page, Nucl. Phys. B443 (1995) 233 and Phys. Rev. D52 (1995) 
1706. 

J. Adomeit et al, Zeit. Phys. C71 (1996) 227. 
J.Z. Bai et al., Phys. Lett. B472 (2000) 200. 



10 



FIGURES 



600 



500 



^ 400 

^ 300 

^ 200 

s 100 



h — ^ — 


1 1 1 1 1 . 




40 




(a) i 










Me 


30 








20 














10 


= 1 ^..UY^ 


1 , = 




n 



0.25 



T — I — I — I — |— 1 — I — I — r-| — I — I — n — I — I — I — I — I — I — r 

(b) 




I . . . . ii 



1.5 



2.5 



ntK^K Kl (GeV/c ) 



0.5 0.75 1 

(GeV/c') 

FIG. 1. (a) the vr+vr^ invariant mass with invariant mass closest to the K'^ mass; (b) KKtt 
mass spectrum. The dark dashed histogram of (b) shows the estimated background in the analysis 
region {KKtt mass below 2.3 GeV). 
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FIG. 2. Dalitz plots for KKtt mass ranges (a) 1360-1560 MeV, (b) 1600-1750 MeV, (c) 
1800-2200 MeV; (d), (e) and (f) show fitted Dalitz plots. 
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FIG. 3. Projections on to (a)-(c) K^K^ mass, (d)-(f) K-k mass for the three mass intervals 
of Fig. 2; histograms show the fit. 
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TABLES 



TABLE I. Branching ratios (BR) of 77(1440) integrated over its width. 



Decay Channel 


BR (%) 


K*K 


0.70 ± 0.05 


kK 


0.13 ±0.03 




0.09 ± 0.03 


pp 


0.08 ± 0.03 


TABLE II. Branching fractions (BF) for production and decay. 


Values are corrected for all 


charge states in KKir. 




Process 


BF (%) 


1) J/ip 7??(1440), ?7(1440) ^ KKir 


(1.66 ±0.10 ±0.58) X 10-3 


2) J/^ 7/i(1285), /i(1285) KK-k 


(0.61 ± 0.04 ± 0.21) X 10 


3) J/^ ^ 7/i(1420), /i(1420) ^ KKti 


(0.68 ±0.04 ±0.24) x 10"^ 


4) J/ip 7?7(1800), ?7(1800) ^ kK 


(0.58 ± 0.03 ± 0.20) X lO'^ 


5) Interference between (1) and (4) 


(0.15 ± 0.01 ± 0.05) X 10-3 


6) Interference between (1) (3) 


(-0.03 ±0.01 ±0.01) X 10-3 


7)J/V' ^ 7??(2040), ?7(2040) ^ kK 


(2.1 ±0.1 ±0.7) X 10-3 
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